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SUPERSONIC  FLOWS  AROUND  A  CIRCULAR  CONE  WITH  OR  WITHOUT  BLOWING  ON 
THE  SURFACE 

E.  Carafoli,  C.  Berbente  and  P.  Marinescu 

SUMMARY 

This  article  studies  the  flow  around  a  straight  circular  cone  with 
or  without  blowing  in  the  presence  of  a  supersonic  current  parallel  to 
the  axis  of  symmetry.  First  we  consider  the  case  of  a  solid  circular 
cone  for  which  a  nearly  exact  solution  is  found  for  the  equation  of 
motion.  Thus  formulae  are  obtained  for  simple  calculation  to  determine 
the  angle  of  the  shock  wave,  the  fields  of  speeds  and  pressures.  The 
results  obtained  are  very  consistent  with  the  data  of  other  authors 
using  numerical  methods  and  with  experimental  data.  This  fact  made  it 
possible  to  extend  the  method  to  the  case  of  flow  around  a  permeable 
porous  cone  with  suction  or  normal  injection  on  the  surface. 

The  results  obtained  may  be  applied  directly  to  the  calculation 
of  the  fuselage  of  supersonic  and  jet  aircraft,  or  the  thermal  pro¬ 
tection  of  the  tip  of  aircraft  by  sonic  injection  and  rockets,  or  the 
thermal  protection  at  the  tip  of  aircraft  by  injection  of  a  liquid  or 
gas.  Likewise  the  use  of  suctions  through  the  walls  reduces  the  inten¬ 
sity  of  the  shock  wave  for  a  certain  period  of  the  flight,  eliminating 
a  series  of  harmful  effects  of  the  latter. 


1.  MOTION  AROUND  A  SOLID  CONE  WITHOUT  BLOWING 

A  circular  cone  is  considered  located  in  a  supersonic  current  of 
velocity  Uoo(fig*  1)  • 


Designating  by  V(0)  the 
of  the  vector,  we  may  write 


velocity  component  according  to  the  radius 
V*  =  V(  •)-!*, 


=  V\ 


d*F 

dO* 


O' 


Since  the  movement 
Taylor-Maccoll  equation 


is  conical,  the  velocity  field  is  given  by  the 

•a 

V't\ 


('  •  (*  -  7,') 


V 


0. 


I'-*) 


in  which  a  represents  the  local  velocity  of  sound,  given  by  the  relation 

fl*  a*  -Y“- (!■'*- (•<) 

where  Y is  the  ratio  of  the  specific  heat  (y=1.4  for  air). 

Introducing  the  velocity  of  interference  V,  given  by  the  equation 


-2- 


P  ^  ~  foa  6,  (v 


Vmt sine,  r"  =  »"  a  i*  (os 


equation  (2)  becomes 


6)  (4) 


0  a*  )V  +  w*  °lg8  +  (2~  “,*)  r  -  « 


Since  it  is  not  simple  to  solve  equations  (2)  or  (5)  we  propose 

2  2 

to  linearize  them  by  replacing  the  nonlinear  term  V*  /a  by  the  function 


t2  (0) 


\r't 

a* 


suitably  chosen,  thus  the  solution  which  results  will  have  a  simpler 
form  leading  at  the  same  time  to  more  exact  results. 

Consequently,  instead  of  equation  (5)  we  will  solve  equation 
(I  -  <*)  e"  |  r'cttfO  (-  —  t*)r  -  0,  (•) 

whose  general  solution  is: 

u~  (<  1  ’  1>]  tos  8’  W 

where  the  coefficients  C  and  D  are  constant  quantities,  while  I  is  the 
result  of  the  following  integrals: 

/  f  d0  ...  ,,  f‘tgO‘10 


f  ™  tW**} 

3 //mm*  o  3  l  - 


From  equation  (8)  we  deduce  immediately 

■—  f  —  tg  0.  (10) 

ii  ids  o  r. 

to  determine  the  constants  C  and  D  we  use  the  boundary  conditions  with 

( C>=9  )  and  for  the  shock  wave  (9=9  )  .  Since  the  speed  normal  to  the 

surface  of  the  cone  is  zero,  we  will  have 

e-ef,  v;  =  r;  -  r.  «iner  -  o.  on 

From  the  equation  of  conservation  of  masses  and  impulses  we  have 
the  relationship  A 

8=0”  r*=0>  M2,4 


(12a) 


r;  =  r.o -t)fi -±' 

|  V 

(12b) 

where 

Jfc  =  ?—  •  K,  =  M„  sin®,. 

(13) 

Y  +1 

Applying  conditions  (11) ,  (12a)  and 

(12b)  we 

then  deduce 

E.co*Q,~{CI‘  +  D)  S,n0(  =  s,n0r» 

(14a) 

Cl,  4®  =  0, 

(14b) 

H.  cos®.  +^)««nG,  =  (l  ~*)(l 

(14c) 

Starting  from  them  we  obtain 


—  =  Hf  -  /.)«»*«  (D  =  —  r/j, 

U  oo 

H{  sin®,  cos6c 
1  -f  //,( /,  —  /, )  sin  0,  cos®, 

si»*0,-  1  + 

il/i  «  (I  -  k)  H, 

V 

* 

r  =('(  1  -  </-  AHii.oV 

C„  \H  OOS0  ) 


(13a) 
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(Hi) 


where  H  ,  H  ,  I  ,  I  represent  values  of  H  and  I  for  0=0  and  0=0  .  To 
c  s  c  s  ^  c  s  ^2 

determine  the  function  t  ,  we  may  observe  that  (fig.  2)  the  function  V  / 

2 

a  and  its  first  derivative  become  equal  to  zero  for  e=ec»  but  the  flow 

after  the  shockwave  has  a  normal  subsonic  component.  Even  in  the  case 

2 

when  the  motion  is  incompressible  we  will  have  t  =0,  and  we  obtain  from 
equation  (9)  //„  sin®,  lu  -  ^  |  In  .  (17) 


We  may  note  the  fact  that  relation  (17)  furnishes  results  close  to 
the  real  ones  which  makes  it  possible  for  us  to  choose  a  simple  expres¬ 


sion  for  t  in  the  form 


<*  =  b(.r-  1), 


..  I  ‘i"e  v.  . 

*  bin  0e  ) 
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(is) 


-4- 


where 


b  = 


t; 


j-.-i 


,,  v?  =  kK'-n-k 
*  (l-f*)tfjf-ifc  ’ 


(19) 


2 

the  values  of  being  determined  on  the  basis  of  the  boundary  condi¬ 
tions  (12a)  and  (12b). 


Expression  (18)  as  resulting  from  the  applications,  gives  nearly 
exact  values  and  leads  to  practically  applicable  expressions.  These 
expressions  are  simple,  being  established  in  the  form  of  a  correction 
to  the  incompressible  solution  (17) . 


Introducing  expression  (18)  in  (9)  we  have 
_  rotgdO  _1_  f  dx _ 

31  -!•"«)  xri- b(x  -  1)1  ' 

The  integral  (20)  is  solved  immediately  and  results  in  the  expres 


(20) 


sion: 


i  =<1-°£T 

If  (sin  Sr* 


(sin  0)®(1  —  nr)' 
sin  6 


(21) 


in  which  some  constants  are  given  up  which  can  be  included  in  the  gen 
eral  constant  C  and  we  note 


To  calculate  the  integral  I  we  observe  that  expression  (21)  may  be 
expanded  in  the  form 

_L  -  H-*1-  sin0)  f(L- _aTl  _  [ 1  si"  6)  J  ( >  -  K  * r)  „ 

H  sin  6  sin  6 

«  — °  -  oS,(1  “  -  i0  -  -®M"  ,  sin-0,  (23) 

sin  9  sin*0e  sin-0,. 

where6o  and  6yare  constants  chosen  suitably  to  represent  better  the  de¬ 
velopment  of  the  two  binomials  of  expression  (23) .  It  is  observed  that 
fo  8„  =  8V  =  0  we  obtain  expression  (17)  and  Hq  .  Consequently  6  and  6^ 
represent  correction  terms  before  the  basic  integrals  and  Ig.  Like¬ 


wise  it  is  possible  to  put  the  function  H  in  a  simpler  form,  retaining 
only  the  important  terms  of  (23)  and  certain  terms  in  l/sin0  and 
terms  in  sin*'1©  which  have  a  larger  coefficient  before  the  final  term. 
Consequently  we  may  write 

J _ 1—8 

u  ~  +  *  »in“_,0. 

«  siD  0  (24) 

where  the  coefficients  6  and  e  are  determined  from  the  condition  that 

the  function  (24)  is  obtained  from  the  points  0=6  and  0=0  ,  where  the 

c  s 

values  of  H  in  these  points  are  given  by  (21) : 

(2r.) 

Hr  V  sin-0,  /  H,  V  sin-0,  / 

From  conditions  0=6  and  0=0  ,  we  have 

c  s 

_L  =  LnS  4-  esih’-'O,,  -L  =,  1  ~  *  f  esin-'1©..  (26) 

H,  sin  0,  %  H,  sin0. 


and  consequently  we  obtain 


H,r,  sin0r  —  //fsin0, . 
H'H.ir,  -  1) 


1 _ 1  -J 

H,  sin--‘0c  sin-0, 


We  emphasize  the  fact  that  the  expression  (24)  is  only  used  to 
calculate  the  integral  I,  the  single  place  where  it  is  necessary  that 
the  deviations  are  the  smallest,  as  follows  from  the  integration. 


With  expression  (24)  obtained  for  1/H  the  integral  (9)  may  be  ob¬ 
tained  easily  in  the  form 


l  ( 1  -  8 )  Iq  -t  e  \ - 

'  cos*0 


(28) 


As  may  be  seen,  the  basic  term  of  the  integral  I  is  Iq,  the  others 
being  correction  terms,  a  fact  which  makes  it  possible  to  simplify  the 
calculations.  Thus  all  the  above  expressions  are  easy  to  manipulate 
and  lead  to  results  which  practically  correspond  to  the  numerical  cal¬ 
culations  . 

The  verifications  carried  out  show  that  the  same  results  are 
valid  for  the  complete  flow  conditions  in  which  the  shock  wave  remains 
attached,  including  the  transonic  conditions,  corresponding  to  n=4  and 
leading  to  an  expression  for  I  deduced  from  (28)  in  the  form 

^  =  (1  —  $)  I0  -f  cfco&O  H - —  ^  .  (29) 

\  leos  0  ) 

Determination  of  the  angle  of  the  shock  wave.  To  determine  the 

angle  0  of  the  shock  wave,  we  will  apply  equation  (15c)  which  leads  to 
s 

the  expression  sjn^  _  1_  _ Hc  sin  0C  eosOr  tgO, _ 

Af*  (1  —k)U,  11  +He(It—Ie) sinft,  cos0, ) 

This  expression  may  be  solved  in  a  simpler  or  more  complicated  man 
ner,  through  numerical  examples.  This  is  accomplished  more  easily  by 
using  the  method  of  successive  iterations  which  leads  but  after  the  se¬ 
cond  approximation  to  exact  values.  For  this  purpose  we  will  seek  the 
most  suitable  starting  expression  for  the  real  value  using  for  H  in  the 
first  approximation  the  expression  given  by  (17),  i.e.  Hc=sin  0  ,  Hg= 

sine  0  .  We  will  write  \  1  +  _ I _ eos^  (31) 

s  Mi  1  —  k  1  +  (/0J  — /or)8iJ)*0c  CO80,  cos 0* 

Starting  from  this  and  observing  that  the  last  two  factors  of  the 
second  term  of  the  right  hand  side  of  equation  (31)  are  mutually  com¬ 
pensated,  tending  to  the  value  1,  we  will  write  in  the  first  approxi¬ 
mation  -iiu*0*  =  ^ 


In  the  above  equations  (31)  and  (32)  the  index*  was  introduced  to 

★ 

designate  the  starting  value  6  After  calculation  of  0*  from  (32)  we 

s  s 

will  introduce  this  value  in  (15b)  and  (21)  and  we  recalculate  the  value 
of  8  from  (15c)  in  the  form 

5 


sin*  0, 


I  f*  Iff0* 
Ml  o  -At//; 


With  this  value  of  9  ,  we  will  calculate  the  new  characteristics 

s 

of  motion  C,  H,  I. 

Calculations  of  pressures:  After  the  pressure  drop  by  the  shock 
wave  (transformation  accompanied  by  losses) ,  the  gas  unde-goes  isen- 
tropic  compression,  reaching  a  maximum  pressure  at  0=8c 

Designating  by  p,  the  immediate  Dressure  after  the  t^ck  wave  we 
have  _  2Y  V*  si«*0,  - Y  (1  .  k)  K}  *•->.,  (3*. 

p,  Y  +  1  *  •  1 

a  ratio  which  we  designate  by  X ,  to  simplify  the  writing  of  the  following 
formulae. 

After  the  shock  wave  the  transformation  of  the  gases  are  isentropic, 
the  ratio  of  the  pressure  p  under  the  current  angle  8  and  the  pressure 
p* may  be  written  in  the  form  Y 

P  _ p  P’  -  ys  “  r’*V  1  MM 


p  p  p,  _  /(IV*  -  Vs  - »  *y 

p.  v.  p »  l « *  -  v;  -  »7 ) 


where  W  represents  the  -^vimum  speed,  given  by  the  equation  /l / 

IV*  I  _  k 

=  1  -f  -  .  (36) 

tl  kUi 


From  the  boundary  conditions  in  (12a)  and  (12b)  we  obtain 

n*  _  v*  -  r:-  i  -  k/.  .  i  -  my* 


n-*  _  v.*  -  17  _  l  -  »/,  t  i  - 
r*  Mi  l  *K;  ) 


In  this  case,  for  a  certain  angle  0  and  for  the  surface  of  the  so¬ 
lid  cr,-">  (-=(4  V=V  ,  V 1  =0 )  we  have 


). 


"  V  ) 

t'V) 


’  i  i Y  1  v.2  (4  -  1 ' )  v  ' 

JL  2  r  n) 

"  1  ‘(‘V)  I 

Knowing  the  ratio  of  pressures,  from  the  equations  (35)  and  (38)  it 
is  possible  to  calculate  the  pressure  coefficients  defined  by  the 
relation 

c,  2  f-e-  -l)-  (40) 

v-Wi  I/**  ) 

Verification  and  comparison.  The  comparisons  are  made  for  two  im¬ 
portant  parameters  of  the  flow,  specifically  the  angle  of  the  shock  wave 
®s  (equation  33)  and  the  pressure  coefficient  on  the  body  Cpc  (equations 
39  and  40)  . 

From  figures  3  and  4  we  observe  that  the  results  obtained  practi¬ 
cally  for  9  and  C  coincide  with  the  values  obtained  by  numerical  in- 
tegration  of  equations  (2) ,  (3) ,  (4) ,  for  a  large  interval  of  mach 

numbers  (r^=l .  25-5 . 7)  and  the  deviations  (8c=10°-50°)  covering  the  en¬ 
tire  supersonic  and  hypersonic  range  are  moderate. 


It  is  observed  that  the  method  is  applicable  also  in  cases  when  the 
flow  conditions  are  transonic,  as  occurs,  for  instance,  at  1^=1. 25, 
9c=20°  or  at  ^=2,  ec  ,=40°.  In  the  latter  case  the  flow  is  very  close 
to  the  point  of  detachment  of  the  shock  wave  (6^  ^^^*41 . 5°)  . 

We  may  therefore  observe  the  satisfactory  consistency  of  the  ana¬ 
lytical  formulae  and  the  results  obtained  by  numerical  integration,  as 
well  as  the  experimental  data,  which  implies  a  very  high  precision  for 
the  distribution  of  velocities. 


2.  FLOW  AROUND  A  POROUS  PERMEABLE  CONE  WITH  SUCTION1 


It  is  assumed  that  through  the  walls  of  a  cone,  with  a  half-angle 
at  the  vertex  6=9p  (fig.  5)  the  suction  is  carried  out  at  normal  speed 
W  •  On  the  vector  radius  9=6C»  the  normal  speed  is  equal  to  zero  and 
therefore  the  cone  forms  with  this  angle  a  straight  definite  shaft  with 
solid  equivalent  by  analogy  to  the  previous  case. 


(41  > 


Designating  by  £  a  coefficient  in  the  form 

ft  i  »»,. 

P  1  ..  •  0  <  0  s:  | 

l  ,>1110,, 

for  the  characteristic  of  suction,  and  taking  the  boundary  conditions 
just  as  in  the  previous  case,  with  the  difference  that  the  walls  are 
porous  we  will  have 


resulting  in 


0  - -0„,  Y'v  -  t\  -  T*  sin 8,  ^  wpt 


0t'xsin  8P 


(4Z) 

(43) 


and  consequently  we  find  for  C  the  expression 

_  p//„  sin  6,  cos  __  H,  sin  8,  cos  8,. _  ^ 

1  4-  //„(/,—  /„)  sin  1  —  /t)  sin  0C  cos 6C 

obtained  in  the  same  way  as  expression  (15b) . 


We  will  consider  three  important  situations  of  the  calculation: 

a.  Direct  problem.  The  angle  9^  and  suctions  Wp  or  the  parameter 
fi(  o<e£i)  are  given,  together  with  the  Mach  number  of  the  flow  M  . 

From  equation  (44)  we  calculate  the  angle  9^  of  the  equivalent  cone, 

o  ^ 

using  for  H  the  expression  H  =sin  9 . 

We  will  obtain  successively,  observing  that  the  denominator  does 
not  differ  much  from  unity, 

1A  similar  article,  "Supersonic  Cone  with  Blowing",  E.  Carafoli  and  C. 
Berbente  was  published  in  The  Aeronautical  Quarterly,  London,  1976. 
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(451* ) 


2siti0®^l  —  siii*0e°  ft  0  sin  0psii)  20„, 
siu*0,°  *  1  -  yf-  0  sin  0pliTu  20p.  (45c  > 

With  the  value  of  9C  given  by  (45c)  the  problem  is  solved  in  the 

case  of  a  solid  cone  and  we  obtain  the  values  of  H®,  H®,  C° ,  1^,  1°. 

pc  sc 

Then  we  calculate  the  expression  for  the  angle  0  in  two  approximations 

C 

of  the  equation  (44)  in  the  form 

sin  2 0,  m  2 +#«(/,»-!«) sin  20« 

sin20„  P  Z/?2  +  fl;(/«-/J)8in20J,  ‘ 

This  new  value  of  0  deduced  from  (46)  is  sufficiently  accurate  and 
no  more  corrections  are  needed.  The  upper  index  (0)  was  introduced  for 
designating  the  first  approximation. 


b. 


Indirect  problem.  The  parameters  B  and  0  are  given,  and  we 

c 


seek  the  angle  0  , 

P 

tion  of  parameter  3. 


which  requires  the  implementation  of  a  definite  suc- 
In  this  case  we  calculate  from  (45b)  the  value  of 


the  angle  9p  through  a  first  approximation,  which  is  designated  by  9p. 
We  will  obtain  just  like  for  (45c) 


8in*0J  — 


1  1 1  —  y  sin  °r  s 


sin  2  0r* 


(47) 


The  correction  of  the  angle  9^  is  accomplished  using  once  again  the 

equation  (44)  in  the  form 

sin  2 0P  =  J_  He  2  +_/§(/,  -  /°)sin20g 

sin  20r 


K  H  HC(I .  -  I r)sin20r 


(4H) 


The  necessary  speed  of  suction  is  calculated  by  means  of  relation 

=  —  1)  sin0p.  (49) 


c.  Mixed  problem.  We  observe  that  the  angle  of  the  solid  equiva¬ 
lent  cone  9  may  be  used  as  a  parameter  for  the  suction.  This  means 


This  means  that  9  may  indicate  the  intensity  of  the  section  by  choosing 

P 

the  angle  0<e  <«  ,  In  this  case  the  flow  is  determined  and  there 
c  p 

only  remains  the  calculation  of  the  suction  speed  for  9=0  from  equations 
(42)  and  (16)  r  i 


r[ffpcos0p  ilp  -Kin6p 


It  is  found  that  it  is  simpler  to  solve  this  in  the  third  case, 


CONE  AROUND  A  POROUS  CONE  WITH  INJECTION 


V'  X* 


V“> 


A  J 


If  through  the  porous  wall  a  fluid  is  injected  at  a  constant  speed 
(Wp>0),  for  example  for  the  purpose  of  thermal  protection,  a  surface 
of  separation  will  occur  between  the  injected  fluid  and  because  of  it 
a  current  to  infinity,  located  under  an  angle  0c> 0p  (fig.  6).  On  this 
surface  of  separation,  the  normal  speed  component  is  cancelled  (V  =0) 
and  may  therefore  be  defined  also  in  the  case  of  suction  as  an  equiv¬ 
alent  cone.  However  we  should  observe  that  the  radial  speed  component, 
that  is  V,  may  present  a  drop  at  9=9c.  This  drop  occurs  even  if  the 
fluid  injected  is  identical  to  the  external  one,  but  it  has  a  different 
stagnation  pressure.  Consequently  the  conical  surface  9=9C  represents 
the  conical  layer  which  maintains  a  stable  form  if  the  pressure  in 
these  two  flows,  internal  and  external  of  the  cone  8„  are  equal  (p  .  = 


The  external  motion  around  the  equivalent  cone  of  half-angle  ec 
has  been  studied  previously  (paragraph  1) .  We  will  now  study  the 
motion  of  the  inside  of  the  equivalent  cone,  which  may  be  treated  inde¬ 
pendently  of  the  external  one.  The  problem  is  therefore  reduced  to  the 
study  of  the  motion  inside  a  solid  cone  of  half-angle  ©c  after  the  in¬ 
jection  of  a  fluid  through  a  porous  cone  of  half-angle  9^  coaxial  with 
the  solid  one  (fig.  7) . 

Flow  inside  a  solid  cone  following  a  subsonic  blowing  through  a 
coaxial  internal  porous  cone. 

Let  us  consider  a  porous  cone  of  half-angle  9=6p  inside  a  solid 

cone  of  half-angle  ©c  (fig.  7) .  The  fluid  injection  is  carried  out 

through  the  porous  wall  along  the  normal,  at  a  constant  speed  w’»  more 

than  0,  along  the  radius.  In  the  following  we  will  assume  that  the 

2 

speed  W  is  subsonic  (t  <  1) .  The  problem  is  bound  with  the  solution 

P  P  2 

of  equation  (2)  which  we  shall  linearize  by  introducing  the  function  t 
in  the  form  (1  ,*>  y  +  F’cteO  +  rj  ..  /*,  y  0.  (.r>D 

The  difference  for  the  previous  case  for  the  solid  cone  consists 

2 

in  the  fact  that  the  function  t  must  be  considered  on  the  descending 

branch  (9<0c  fig-  2)  and  the  speed  of  sound  is  obtained  from  the  expres 

sion  v  _  1 

a*  -  a2,  +  y (  V-  \  H  *  -  V*  -  I"*),  (®*> 

2 

where  9p  represents  the  speed  of  sound  0=9^. 

The  general  solution  of  equation  (51)  is  similar  to  that  given  for 
equation  (8)  .  We  will  write  .. 

H  l  =  [01  +  F)  cos  e,  (.53) 


where  the  function  I  is  given  by  (9) ,  by  the  constants  G  and  F  are  de¬ 
termined  from  the  following  boundary  conditions 


We  obtain  immediately  from  these  two  equations  (54) ; 


H'H,  Mn20j_*ii«a0„ _ 


,  _ . _ _ _  .  .  —  •  (.">r».») 

6,  2Hr  sin  20,-2 H,  *in  20r-f  (/,—  f„) 11  20(-s,n  28, 

F  =  0f - 2 - /rV  (•’»•»>> 

\  H,  «in  2  0,  / 


For  the  function  t  we  will  choose  an  expression  of  the  form 

r  =  ( *in0  .  m  -1,2,3 .  (r>6> 


-OJ 


which  is  similar  to  the  form  (18)  but  for  negative  n. 

The  coefficient  b  is  determined  from  the  equations 
h  -  J*:vp  j s  _  * H'J 


b  = 

l— t„ 


For  t  =0,  we  obtain  expression  (17)  for  the  functions  HQ  and  IQ 


the  form 

H0  =  sin  0,  I0  =  1  -f  In  tg  •  (58) 

cos  0  2 

For  this  case  which  corresponds  to  an  incompressible  motion,  the 
constants  which  we  shall  designate  as  FQ  and  GQ  are  obtained  from  the 
expressions  (55) .  The  same  movement  may  take  place  if  a  liquid  is  in¬ 
jected  for  cooling  purposes.  Let  us  consider  the  angle  8=9  corres- 

2  1 
ponding  to  t^=l  (fig.  2)  where  6^  represents  the  boundary  value  which 

is  never  reached  in  subsonic  blowing,  that  is  we  will  always  have 

9^.  Taking 

a  =  .  mv  =  1  -  «,  (r>9) 


/sin  0,\"  b  , 

— - *  I  =  -  -  ,  wv  =  l  - 

\  gin  0,/  1  +  b 


we  can  write 


J  —or  d.r 

>’1  J  JT  — 


f  dr  const. 

Jt-(j  '  // 


or  choosing  suitably  the  constants 
l  1  1 


1 

1 

[  1  — « 1  —  sin  0)  | 

ft. 

sin"  0, 

V  Rin0 

/  s.n-o,  \ 

l 

Kili"  0 

1 

\  sin-0  ) 

-15- 


Since  v  and  o  have  low  values,  the  expression  (61)  may  be  expanded 
in  suitable  manner,  which,  just  like  expression  (23)  will  be  used  only 
for  integration. 


Thus  using  suitably  chosen  constants  ,  6^ ,  we  may  write 

I  1  1  -  8,(1^  sii.O)  1  -  8 B  , 

II  sin 6  ]_§  Hiu*  ®i  ^ sin 0  ° 

sin"  6 

j.  -8<8osin"01  +  28.(1-  8„)  sin"  0, 

sin-tt  (siu  0)"’ 1  '  1 

Just  as  in  the  previous  case  we  can  consider  a  suitable  value  for 
6oand  6^  for  the  purpose  of  a  better  representation  of  expression  (61) 
Observing  that  the  first  and  last  terms  of  (62)  are  preponderant,  the 
above  expression  may  be  written  with  sufficiently  good  approximation  in 
a  form  similar  to  (24):  1  l  —  8' 

H  sin  0  (Sin  0)«+i  * 


The  constants  6 ’  and  e 1  are  determined  from  the  conditions  that  at 

the  ends  of  the  interval  e=0p  and  e=ec»  the  expression  (63)  coincides 

with  (61)  that  is  H  will  have  two  values  H  and  H  . 

P  ^  •  m  w 

//,  (sin  0cm  -  nv  (Sin  O,,)"**  f  i  _  slU"0'  V.  (64) 

\  sin"0„/ 


By  this  means  we  may  easily  deduce  that 


j  g,  //,,  >in  0,  -  H.j'f  sin  0, 
H  H  .i  l  -  x„) 


i"e{Tf_(l“s')]  (tt5) 


Consequently,  the  integral  I  will  have  the  form 


/  =  (1  -«')/„  +  e'(  - 

J  tos* 


(10 

0  sin"* 1  0 


Results  very  near  to  exact  ones  are  obtained  for  m=3  when  the 
expression  (66)  leads  to  ,  /  1  \ 


Calculation  of  pressures.  We  designate  that  the  pressure  on  the 
internal  cone  0=9^  and  pc  the  pressure  on  the  external  cone  6=6C.  Desig 
nating  by  p  the  pressure  at  a  current  point,  we  will  have 

t  »?- ti.-r.'.-n-1.  * y!*.  m 

p,  ur>  -  n  -  ;  p,  \  IV.-  V;  v,4 


where  the  maximum  speed  W  is  given  by  the  expression. 

i 

wf  =  \’i  4  r;*  +  (69) 

Y  ~  1 

From  here  we  calculate  the  pressure  coefficient 

C„^  2  (V  -\).  ('«>) 

Y Ml  \pP  I 

M  being  the  Mach  number  at  9 =0  . 

P  P 


Determination  of  the  lines  of  current.  It  is  useful  to  know  the 
form  of  the  lines  of  current  in  the  region  inside  the  two  cones  in  view 
of  multiple  applications.  In  the  first  row  we  will  give  the  expression 
of  the  angle  a  of  inclination  of  the  total  speed,  on  the  side  of  the 
surface  of  the  porous  cone 

tp*  =  ?>_  -  W *i"  2%)-2l(Hr  sin  20r)  4 
V,  2I(H,  sin2er)-(fe— /,) 


We  will  determine  below  the  function  of  the  current  V  using  the 
equation  of  the  lines  of  current  in  spherical  coordinates,  for  the 
special  case  of  conical  movement 

d  R  Rd0' 

=  - • 

V  V 


The  function  of  the  current  is  introduced  in  the  form 

$  =  R*C  -£  sin  erd0  =  R*  (*  (01  4-  F)  *inO  d0,  (73) 

V  Pl  ’*c  Pi 


where  for  e=ec  we  obtain  the  line  of  current  ^=0,  while  represents 
the  reference  density  and  it  may  be  shown  that  p.=p  .  In  case  of  in- 


compressible  flow,  p=p  ^constant ,  we  have  t  =0  and  we  obtain  for  the 
function  of  the  current  the  following  expression 

J?*£  ((V.4-F,)  sin  8d8  =  -^<70|sin*  8  In  tg  ^  -  cos  ej  4 

4-  Fo«in*0  4-  const  j*  1*4) 

where  the  constant  is  determined  from  the  condition 

For  the  general  case  on  the  basis  of  the  expressions  (67)  and  (73) 
observing  that  we  can  find  for  the  average  value  and  therefore  the 
ratio  p-'p^,  of  small  variations  (p/p^*  1),  we  have 

4,  =  —  —  J')  sin*01n  tg^  -  cos  e| 

+  —  (7e'(tg0  —  2sin 20)  4-  F  sin‘0  ‘  const. |  ("•'*) 

Thus  we  obtain  the  analytical  expression  for  the  function  of  the 
current,  which  may  be  applied  without  any  difficulty. 

Examples  of  calculation  and  choice  of  the  injection  parameters. 

In  case  when  we  consider  a  motion  in  two  coaxial  cones  (Fig.  7) 
of  which  we  know  the  half-angle  of  th'  solid  ^ud  0=  0,.  and  the  half¬ 
angle  of  the  internal  permeable  cent  0  =  0P  .  . etuti  wilh  the  speed  of 
sound  av  te  level  of  the  porous  wa.:  —  8pi  ioi  the  blowing  speed 

W  alonK  the  normal  X?-  may  calculated  along  the  walls  according 
p  6  ’  P 

to  the  relation 

I*  -  *-?  <«*<!>. 

0; 

which  is  implemented  if  the  blowing  is  subsonic.  In  this  case  all  the 
relationships  of  the  calculation  are  explicit  and  there  is  no  difficul¬ 
ty  in  their  application. 

In  case  the  fluid  is  injected  inside  a  current  from  infinity 
(fig.  6)  the  cone  8=9c  represents  the  surface  of  separation  between 
two  flows  (one  internal)  (9^  9^)  and  the  other  external  (9>  8  ) .  We 
may  disregard  the  eddy  effect,  located  at  e=e  ,  as  well  as  the  heat 


exchange,  because  the  latter  is  generally  of  secondary  nature  and  rep¬ 
resents  only  a  correction  which  may  be  added  subsequently.  The  prob¬ 
lem  is  solved  by  considering  simultaneously  two  flows,  the  external 
one  whose  solution  has  the  form  (8)  and  the  internal  one  given  by 
equation  (53) . 

The  stabilization  of  the  quantities  characterizing  the  flow  leads 
to  the  following  situations  of  calculation: 


a.  Direct  problem.  The  characteristics  of  the  external  current 

are  given:  U<*  M«,  p  ,  angle  9  ,  the  speed  of  sound  a  ,  the  speed  of 

»  P  p 

injection  Wp  and  pressure  Pp,  and  we  seek  how  to  determine  the  position 

of  the  equivalent  cone  ec  and  the  position  of  the  shock  wave  0  ,  at 
the  same  time  as  the  entire  field  of  motion.  Together  with  the  con¬ 
stants  C,  D,  G,  F  we  have  six  unknowns  for  which  we  have  six  equations 
given  by  the  two  conditions  9=9  (equations  12a  and  12b) ,  three  con- 

o 

ditions  for  9=9c  (equation  11  and  54b  for  the  external  flow  and  the 
internal  flow  and  equality  of  pressures)  and  the  condition  9=9p. 

Taking  into  account  the  expressions  of  the  integration  constants  (15a, 
15b,  55a  and  55b)  the  problem  is  reduced  to  finding  the  angle  of  the 
solid  cone  9  and  the  angle  of  the  shock  wave  9  from  the  system  of 
equations : 

sin2  6.=  *  -f  -  ...  , - //„ sin  2 6C _ 

Ml  <!-*,)//„  2  +  H „(/„-/„)  sin  20,  ’  '  ’ 


sin2  6,  =  1 

+  .  Jg0. 

Ml 

<»-*,) /i, 

i  +*-1Jfi 

(i- 

2 

/,  i 

—  if  \ 

i  I  k,  -j-  - 

.1 

l 

*.*  J 

where  the  indices  e  and  i  refer  to  the  external  and  internal  flows 
respectively.  From  expression  (77)  we  observe  the  appearance  of  a  new 
important  parameter;  the  ratio  p  /p  . 


The  system  of  equations  (76)  and  (77)  is  easily  solved  by  attempts 

observing  that  at  the  value  chosen  for  the  angle  6  (6c  6  ),  both  the 

c  p 

internal  and  external  flows  are  completely  determined.  This  initial 
value  is  modified  until  equation  (77)  is  verified. 

b.  Indirect  problem.  The  quantities  are  given  which  characterize 

the  internal  flow  a  ,  W  ,  p  ,  9  and  the  angle  8  and  we  attempt  to 

pp  pp  c 

determine  the  characteristics  of  the  external  flow  with  surface  of 
separation  at  8=8c. 

Since  the  internal  flow  is  totally  determined  and  therefore  the 

angle  6c  is  given,  there  remains  the  determination  of  the  constants  C 

and  D  along  with  the  parameters  of  the  external  flow  Utr>,  pw,  and 

the  angle  8  ,  that  is  six  unknowns.  In  this  case  we  have  four  equations 
s 

given  by  the  two  conditions  at  8=0S  (equations  12a  and  12b)  and  two 

conditions  at  8--8c  (equations  11  and  77)  .  The  result  is  that  two  of 

the  parameters  of  the  external  flow  must  be  chosen.  If  we  give  M 

and  U  then  the  angle  8  is  calculated  from  the  equations  (15c)  or  (76) 

s 

while  the  pressure  p  and  therefore  the  stagnation  pressure  from  the 
external  flow  are  obtained  from  (77)  (before  and  after  the  shock  wave) . 
In  this  case  in  the  indirect  problem  the  equation  of  calculation 
becomes  explicit.  We  may  observe  'that  the  speed  U  may  now  be  deter¬ 
mined  from  the  equations  referring  to  the  external  flow,  as  represent¬ 
ing  only  the  stages  of  movement.  Consequently  in  the  indirect  problem 
we  cannot  consider  the  values  and  pm  as  known  simultaneously.  If 
the  values  are  given  for  p*  and  \Jn ,  then  Mw  and  the  angle  8  are  deter- 
mined  from  trials  by  means  of  the  system  of  equations  (76)  and  (77) . 

Comparison  with  other  methods.  Figures  8  and  9  show  the  diagrams 

of  the  results  obtained  on  the  basis  of  the  analytical  methods  proposed 

and  are  compared  with  values  determined  by  numerical  integration  by 

Emanuel  /2/.  We  find  a  very  good  consistency  both  as  regards  the  angle 

of  blowing  a, (formula  71)  and  the  ratio  of  pressures  p  /p  .  The  para- 

2  P 

meters  chosen  were  8  ,  t  ,  and  8  for  which  we  considered  three  values: 


—  Teoria  prtrentd  (  Si ) 
O  Mttodo  numerica  [*) 
(b) 


20' 30'  40'  SO'  Q>60' 


Figure  8:  Key:  (a)  present  theory;  (b)  numerical  method  /2/, 


/ P  •  —  Teorio  prezenld 

U  1/  - — 1  O  Metodo  numerica  [z] 

,  ,  (b) 

l°J£.  jgi  JO*  40'  SO'  e*  60* 

Figure  9:  Key:  (a)  present  theory;  (b)  numerical  method  /2/. 

For  these  values  values  of  a  were  found  according  to  (71)  and  p  / 

Pc  (formula  77)  while  y  =y.=1.4.  The  consistency  is  good  also  in  the 

e 

case  of  larger  values  for  t  . 

P 
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